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oo ' Abstract 
m 

We propose a model for cell migration where epithelial cells are able to detect trajectories of other 

(N 
O 



cells and try to follow them. As cells move along in 2D cell culture, they mark their paths by 
loosing tiny parts of cytoplasm. Any cell moving on a surface where other cells have moved before 
faces a network of cell trajectories, which it tries to restrict its motion onto. With the Tsallis 



modification of classical thermodynamics one can solve the relevant Fokker-Planck like equation 

o 

and obtain experimentally testable distribution functions. We compare the model to experimental 
data of normal mammary epithelial cells and cells which have been genetically manipulated to 
$_i ■ change their degree of cell-cell interaction. 
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INTRODUCTION 



A generalization of classical Boltzmann thermodynamics as recently introduced by Tsallis 
I], |2|, |3| has offered a beautiful way of describing the statistics of a multitude of physical 
systems, which generally involve some sort of long-range interactions, ranging from astro- 
physics j|, |5|] to low dimensional dissipative chaotic systems 0. Systems described by this 
thermodynamics are often referred to as " non-extensive" . This same modification of classi- 
cal thermodynamics also seems to offer a natural understanding of the subject of anomalous 
diffusion. Anomalous diffusion can be thought of the propagation of particles which are 
restricted to move on a certain subspace of full two or three space, such as for example the 
propagation of water molecules in a sponge, or the flow of crude oil in porous soil. In this 
work we argue that the motion of living and interacting epithelial cells can be seen as a 
variant of anomalous diffusion on a network which changes over time. We support the idea 
by experimental evidence. Our experimental findings moreover strongly support very recent 
work on movements of interacting Hydra cells @. 

A living epithelial cell e.g. in two-dimensional cell culture is not a static object but moves 
continuously, a process which is called cellular migration || [J. Cell migration is of vital 
importance in physiological processes such as wound healing or understanding cancerogen- 
esis. For example, intact mammary gland tissue is stable during phases of milk secretion. 
In residual phases of female estrous cycle, however, it is subject to structural remodeling 
where cells become dynamic entities. They proliferate and migrate, eventually forming new 
functional secretory units [JTDj, [II], [T2|] . Dysregulation of these events is believed to contribute 



to cancerogenesis [13, 14, [Lq]. In-vitro experiments preserve important aspects of this ep- 



ithelial biology |lq, 17 , and migration of single cells has been studied in cell culture systems 



on molecular and mechanical levels 1|, 18 



Single cells are by no means passive movers and governed by chance, but influence and 
control their movements themselves. For these movements to occur morphological changes 
are necessary which are realized in the formation of distinct subcellular domains. In Fig. la 
a single cell is shown during movement. After the formation of a so-called lamellipodium, 
(a V-shaped protrusion) this part of the cell attaches to the cell culture surface. The cell 
body is then pulled towards the front part ("leading edge") via the activation of mechani- 
cal intracellular linkages. Eventually the "trailing" edge, a thin, elongated structure at the 
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opposite side of the cell with respect to the lamellipodium, detaches from the surface which 
allows a net translocation of the cell. In this "polarized" fashion the cell apparently crawls 
on the surface of the cell culture substrate, continuously readjusting its direction. During 
movement cells loose tiny parts of their cytoplasm, marking at least parts of their trajec- 
tories. It is known that cells change their movement patterns when the activity of proteins 
involved in cell adhesion is regulated up or down ward. When these proteins for example 
are suppressed, cells loose their ability to migrate in a polarized fashion [19|, |20|. As several 
recent works on velocity distributions of cells show, it seems very likely that cell movement 
depends strongly on the surrounding of individual cells. Experimental results range from 
normal diffusion [21] to exponential velocity distributions [p2 |. It is the aim of this paper 
to show that an anomalous diffusion approach, already discussed in |7]], might be general 
enough to explain the different experimental results on a common theoretical ground. 



THEORY 



To describe the dynamics of a single cell which has the ability to detect the presence of 
trajectories of other cells, we propose the following model. Imagine each cells paints its trail 
onto the surface. The paint slowly fades, and can not be seen after a certain time. Another 
cell which detects the line of paint tries to follow it as good as possible. From time to time 
the cell will be distracted from strict path following by random ("thermal") events. At the 
same time this cell also marks its trajectory with paint. As many cells move across the 
same surface there will evolve a network of paths. Since paint fades and new trajectories are 
added to the network, the network will undergo slow dynamical changes. Any cell equipped 
with a trail following mode will try to restrict its movements onto the network of trails 
of former cells. The situation here is much alike diffusion in porous media where particle 
trajectories are restricted onto certain pathways, and can not propagate through full two- 
or three-dimensional space. In the case of cell migration an individual cell tries to follow 
the dynamical "web" of trajectories other cells have taken before. 

Modeling cell migration as a stochastic diffusion process of this kind clearly involves 
non-linearities in the associated Fokker-Planck like equations 

Ft p{x ' t] = ~iL {F [p] (x ' t)p(x ' t))+ \i^ (jD [p] (x ' t)p{x > t)} • (i) 
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To incorporate the idea that cells follow the gradient of metabolites emitted by other cells 

along their trajectories, a non-linear flow-term of the form F \p) = -^p would be reasonable. 

However, we neglect the flow-term in the following, since a posteriori we find almost no 

experimental evidence for it. To incorporate the anomalous diffusion idea, we decide the 

functional dependence of the diffusion parameter to be D[p]p = p u . With this choice equation 

([!]) can be solved by following the entropy approach by Tsallis |T], [|. Here the the classical 

entropy definition S — — j p(u) lnp(u)du is modified to 

_ 1- Jp(u) q du 
bq ~ q~^l • {2} 

By maximizing S q while keeping energy fixed it can be shown that the resulting p(u), i.e. 

p q [x, t) - ^ [6) 

also solves Eq. ([!]) without the flow-term ^3f. Here Z q (t) is the generalized q partition 
function. The relation of the Tsallis entropy factor q, and v in Eq. (jl]) is given by q = 2 — v 
0. Note that this distribution is a power law, and only in the limit q — ► 1 the classical 
Gaussian pi(x,t) = exp ^^[(t) ~ ~ ^ s recovere d. It is relatively straight forward to derive 
the corresponding velocity distributions of the associated particles 

(3mv 2 



p{v) = B, 



i-(i- q y 



1-8 

(4) 



which has been derived for the first time in p4| , and which reduces to the usual Maxwell- 
Boltzmann result p(v) = (^r) 2 ex P ~ m ^ ne Q — * 1- Here ^ can be seen as an 
inverse mobility factor, meaning that large values of are associated with small average 
cell velocities. If cell migration is in fact describable as anomalous diffusion governed by 
a Fokker-Planck equation like Eq. (]T|) with the special choice of D [p], the cell velocity 
distribution function should follow the form of Eq. (Pi). 

To further characterize the nature of the diffusion process we will also look at normalized 
temporal autocorrelation functions AC(r) = ^|y(C(^)C(^ — r ))*) an d at two-point correlation 
functions C(r) = ([£(£ + r) — £(t)] 2 )t, of an underlying process £(£), where (.) f denotes the 
average over time t. For (fractional) Brownian motion, the correlator scales as C(r) ~ t 2H , 
where H is sometimes called the Hurst exponent. For H = 1/2 the classical diffusion result 
is obtained (Brownian motion), H > 1/2 indicates correlated diffusion, which we expect in 
our case. 



4 



EXPERIMENTAL DATA 



To record the trajectory of an individual mouse mammary gland epithelial cell over 



long time periods of up to 72 hours, epithelial cells [25] were seeded sparsely onto a 2D 
plastic culture dish. The dish was then placed into a mini-chamber mounted on an inverted 
phase contrast video-microscope. This provided incubator conditions and periodic image 
acquisition in 90 seconds intervals was performed f2"Ef . About 900 to 1000 (video) frames 
were taken per measurement. For our purposes we only take the first 25 hours of the 
recording since after that cell density becomes high, and motion of individual cells restricted, 
see video ( Jhttp: / / is2.kph.tuwien.ac.at/thurner / (JFLLMU VF| ) . As seen in the videos, cells 
move within limited regions for many timesteps, crossing and following its own and other 
cell's trajectories. 

From the video sequences the motion of individual cells was extracted semi-automatically 
by a commercial motion tracking routine ^7j. Cell nuclei were tracked based on a convo- 
lution kernel of a 100 square pixel field (about half the size of the cell nucleus at the given 
resolution). These extracted cell trajectories, see lines in Fig. lb and c, were recorded as 
the x t and y t positions in two dimensions (the subscript indicating discrete time). 

Due to the limited data length it is more useful to study velocity distributions rather than 
spatial distributions directly. We work in polar coordinates and transform the trajectories 



{x t ,yt) to a velocity (length-increment v x (t) = x t — x t -i) process v t = ^Jv x (t) 2 + v y (t) 2 and 
a direction process a t = arctan(xt — x t -i, yt — yt-i) €= [— 7r, tt]. 

As a control for the plausibility of our model we genetically modified cells from the same 
cell line in two virus experiments [^, |29|, |30|. The basic idea here is to alter the cell's 
motion program such that its interaction with other cells and hence its ability to follow 
other cell trajectories becomes effectively reduced. In the case of non-interacting cells one 
should expect the model loose its validity, i.e. that cell movement would approach classical 
diffusion, i.e. q — > 1. In the first experiment, we infected cells with a virus carrying 
the coding sequence of enhanced green fluorescent protein (GFP). Since viral infection has 
been reported to significantly influence morphology and molecular homeostasis of mammary 
epithelial cells |32], |33|, , we expected alterations on cellular migration patterns as well. 
This condition will be referred to as " infected" in the following. In a second experiment we 
inserted the coding sequence of gene TPA-Inducible Sequence 7 (TIS7) instead of GFP into 
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the cells. TIS7, a co-repressor of gene transcription and functionally associated with loss of 
epithelial characteristics like cell-to-cell adhesion was chosen to additionally test for eventual 
alterations in cell mobility p9[ . This experimental condition will be called "TIS7" in the 
remainder of the paper and contains consequences of the viral infection itself and the effects 
of over-expressed TIS7. 

RESULTS 

For a first demonstration of correlated motion to be present, we compute C(r) for the 
individual x t and y t components of cell trajectories. We find a Hurst exponent of about 
H ~ 0.7 ± 0.02 for all the experimental conditions (normal, infected, and TIS7), indicating 
correlated motion, see Tab. 1. 

We find the direction process a t practically equally distributed over the interval [— tt, tt]. 
By looking at normalized autocorrelations AC of at, significant (99% significance level) 
coefficients are found which decay as a power law. The situation is shown in Fig. ^|a for 
all three experimental conditions. Power decay exponents were fit by a least square routine 
in the region indicated by the straight lines; the corresponding values are collected in Tab. 
1. Errors (standard mean errors over 5 independently moving cells) in the fit-regions are 
somewhat larger than symbol-size and were omitted for clarity of the plot. Normal cells 
show lower correlations, which decay slightly faster than those of the manipulated cells. 

The normalized autocorrelations for the velocity process v t are given in Fig. 0b. Cor- 
relation coefficients are statistically significant and in the same order of magnitude as the 
direction series'. Again, the tails of the correlations follow a power law; for the correspond- 
ing exponents, see Tab. 1. Here normal cells have higher correlations than infected ones. 
Infected and TIS7 cells show a slower decay than normals. Errors here are about twice 
symbol size. 

In Fig. [5] we show the velocity distribution (histogram) of the normal cells in a log-log 
plot. The solid line is a three parameter fit to Eq. |J fit-parameters for q and for all 
conditions are listed in Tab. 1. At the low limit of v the distribution approaches a constant, 
while the right tail clearly decays in a power law fashion ~ v ~ 7 . We find infected and TIS7 
cells to move faster on average, as seen in the mobility factor or directly in (v t ), see Tab. 1. 
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DISCUSSION 



We provide clear evidence that cell migration of interacting cells is a correlated diffu- 
sion process. This finding is supported by several facts found in our experimental data. 
Correlated motion is shown by the presence of a Hurst exponent significantly larger than 
0.5 and by non- vanishing autocorrelation functions in the velocity and direction processes. 
This finding is in nice agreement with a recent and similar experiment J7|, where an H of 
0.62 ±0.05 was found. The direction autocorrelation functions are found to decay in a power 
like manner, with exponents consistently smaller than 1, opening the possibility for super- 
diffusion. This obviously excludes naive autoregressive or moving average models (ARMA) 
to be realistic candidates for interacting cell migration since these models have exponential 
decays of autocorrelation functions. This finding means that the direction of cells is a "long 
memory" process. For velocity autocorrelation functions the power decay is smaller but still 
present for the normal cells. The value of the correlation coefficient at r = 1 is about the 
same size as for the at process. 

The virus effect is seen in the following way. The most obvious difference of normal and 
manipulated cells is seen in their mobility, which is drastically enhanced for manipulated 
cells. This is expected because we believe that there is less effective adhesion proteins present 
in the infected and especially in the TIS7 cells. Very similar results have been reported in 
0, however, with a different experimental strategy: There cells propagated in different 
aggregates of cell types with different adhesion properties to Hydra cells. 

Infected and TIS7 cells show larger correlations in the direction process, meaning that 
directions are changed less frequently than normal cells do. This is expected since these 
cells detect trails less well and lose their ability to follow other cells. Once moving they 
keep their direction more strictly than normal cells. For the velocity process we observe 
the opposite effect: normal cells have larger autocorrelations, i.e. normal cells change their 
velocities less frequently than infected and TIS7 ones. We attribute this effect to a more 
controlled movement of normal cells when following a path, or equally, a stronger influence 
of random effects on infected or TIS7 cells with damaged movement control. 

We now turn to the main result of this paper. The velocity distribution follows the 
general form of the anomalous diffusion model. Maybe most importantly this model allows 
to naturally understand the presence of the power tail to the right of the distribution. The 
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power can be related (via q) to the v exponent in Eq. (|l|). For normal cells we find q ~ 1.35 
which goes down to q ~ 1.20 for the infected, and to q ~ 1.11 for the TIS7 cells. In the 
Tsallis thermodynamics framework this means, that infected cells approach the classical 
thermodynamic regime, while normal cells with their complex motion control system are 
governed by non-extensive statistics. In this framework it is clear that for non-interacting 
cells ordinary diffusion and Gaussian or exponential velocity distributions were observed in 
PH, |22| , being extreme cases of the present model. If we accept that in our experiments 
infection diminishes the degree of interaction between cells and their trajectories, we see 
that the q value should serve as a sensible quantifier of the degree of interaction. Again, in 
|7j very similar behavior of 2D moving patterns of Hydra cells were found. They reported 
q ~ 1.5, or a power decay in the velocity distributions of about 7 ~ 4 and suggest a different 
way of explaining the movement patterns based on granular flow models. Furthermore, 
these and our findings clearly exclude a Levy flight interpretation of cell movements, and 
led by the observation of trail-following tendencies of cells we believe that cell migration is 
governed by a realization of anomalous diffusion. 

In conclusion, we think that the presented method provides a useful tool to qualitatively 
understand different movement patterns of realistic cell migration, and to quantify the de- 
gree of interaction by the value of the Tsallis q. Since epithelial cell migration plays such 
an important role in situations like wound healing of skin epidermis or cancer metastasis we 
speculate that the above way of quantification could even be a practical pathological char- 
acterization of cells, regarding diverse biological dignities and pharmacological responses, in 
addition to molecular and mechanical parameters. Finally, we note that the basic idea be- 
hind the presented model could also be applicable to other fields, like electrical flow through 
networks of biological neurons. There the effective network topology is also subject to con- 
stant changes due to synaptic plasticity (Hebb or other type of learning). Electrical current 
will flow on excitable media restricted to a given network and by doing so, will re-shape it 
by changing synaptic weights. Here we would predict power law distributions of electrical 
current changes, in analogy to the v distributions presented here. 
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by the Austrian Science Foundation grant FWF, P13577-GEN. 
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Measure 


normal cells 


infected 


TIS7 


Q 


1.350 


1.201 


1.110 


~T~ 


1.20 


0.70 


0.30 


(vt) 


38.8 fim/h 


51.7 fim/h 


71.7 fim/h 


H 


0.686(24) 


0.700(18) 


0.739(31) 


AC decay exponent: v t 


-0.095 


0.026 


-0.009 


AC decay exponent: a t 


-0.713 


-0.701 


-0.591 



TABLE I: Values for the Tsallis q, mobility factors ^r, and average cell velocities (v t ), for the three 
cell types. Data was obtained by LS fits to the aggregated velocity distributions. Additionally, we 
show mean and standard error (from averages over individual trajectories) for the Hurst exponent 
and fits to power decay exponents of the tails in the autocorrelation functions of the cell velocities 
vt and the direction process at- 
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FIG. 1: Normal mouse mammary gland epithelial cells migrating on a two-dimensional culture 
surface. Large magnification of a single cell (a). The leading edge can be recognized as a fan- 
shaped lamellipodium (small arrows) at the bottom of the panel. The trailing edge is marked by a 
star. The direction of migration is indicated by the fat arrow, (b) Extracted trajectory of a single 
normal cell over 1000 timesteps (25 h). 

FIG. 2: Autocorrelation function of the direction process at of cell migration for normal, infected, 
and TIS7 cells (a). Infected and TIS7 cells have somewhat higher values in the correlations and 
show a softer decay, (b) Autocorrelations for the velocity process vt- For infected and TIS7 cells 
the coefficients are lower and the curves are nearly constant, while normal cells again show a power 
decay. 

FIG. 3: Cell velocity distribution (histogram) for normal cells. For small v the distribution 
approaches a constant and decays for larger v in a clear power law. The line represents a fit to the 
Tsallis velocity distribution. The virus experiments are not shown for clarity; their fit parameters 
are gathered, however, in Tab. 1. 
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